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Abstract 

Using a sample of 3.7 x 10^ T^iS) BB events collected with the BABAR detector at the PEP- 
II storage ring, we search for the electroweak penguin decays B"*" — > K~^e'^e~ , B"*" — > K~^fi'^iJ,~, 
B^ — > K*^ e~^e~ , and B" — > K*^ We observe no significant signals for these modes and set 

preliminary 90% C.L. upper limits of 
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1 Introduction 



The rare decays B K£~^i~ and B — > K*l^i~, where i is either an electron or muon, are highly 
suppressed in the Standard Model and are expected to occur via electroweak penguin processes. The 
internal loops in such processes lead to a 6 — > s transition via an effective flavor-changing neutral 
current interaction. Standard Model predictions for these decays are dependent to some extent 
on the modeling of the strong interactions in the hadronic B — > K^*^ transition, but calculations 
indicate that B{B Ki+£-) ^6x 10"^, while B{B K*l+£-) 2 x IQ-^ [§. These processes 
provide a possible window into physics beyond the Standard Model, since new, heavy particles 
such as those predicted by SUSY models can also enter the loops, resulting in significant changes 
to both the decay rates and kinematic distributions 

Experimentally, the small expected rates make searches for these modes difficult. In addition, 
backgrounds from tree-level processes, such as S ^ J/tljK^*\ with J/ip ^ £~^i~, complicate the 
analyses. Searches from CDF Q and CLEO Q have so far yielded only upper limits, although 
in the case of i3 — > K* the limit is only about a factor of two above the Standard Model 
prediction. 

In this paper, we report the results of a preliminary analysis to investigate the backgrounds and 
the ability of the BABAR detector to reject them. We have analyzed an on-resonance data sample 
of 3.2 fb~^, representing about a third of the current BABAR T(4S') integrated luminosity. 

As in many other searches for very rare decays, there are important issues that we face in 
designing event selection criteria in a manner that does not bias the determination of the signal 
yield, or, at worst, create a false apparent signal in a small statistics event sample. The main 
goal of our study is to test the performance of a "blind" analysis in which the event selection is 
optimized without use of the signal or sideband regions in the data. The events in these regions 
are therefore selected without bias, and the sidebands can be used for background estimation. The 
decays B"*" — > J/ipK'^ and B^ — > J/ipK*^, which have the same topology as the signal, provide 
a useful control sample to compare data with Monte Carlo simulation. These events can be used 
to validate the efficiencies predicted by the full Geant321 Monte Carlo simulation. An additional 
control sample is provided by events with K^*^e^fj,'^ candidates. Such events monitor the level of 
combinatorial backgrounds present in the sample. This analysis provides a basis for extending the 
measurement to the much larger data samples that BABAR will obtain in the future. The 3.2 fb^^ 
data sample is comparable in size, however, to that used in an earlier CLEO analysis, and we are 
able to obtain results of comparable sensitivity for the B"*" K~^i^i~ modes. 

We analyze four decay modes: B+ K~^e~^e~, B+ K^fi^/j,^, B^ — > K*^e~^e~, and B^ — > 
K*^ fi~ . The K*^ is reconstructed in the K^ir^ final state. In each case, we include the charge 
conjugate mode as well. 

2 The BaBAR detector and data sample 

The BABAR detector is described in detail elsewhere . All components of the detector are used 
for this study. Of particular importance for this analysis are the charged tracking system and the 
detectors used for particle identification. At radii between about 3 cm and 14 cm, charged tracks 
are measured with high precision in a five-layer silicon vertex tracker (SVT) . This device comprises 
52 modules built from double-sided AC-coupled silicon strip detectors. Tracking beyond the SVT 
is provided by the 40 layer drift chamber, which extends from 23.6 cm to 80.9 cm. Just outside 
the drift chamber is the DIRC, which is a Cherenkov ring- imaging particle identification system. 
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Cherenkov light is produced by charged tracks as they pass through an array of 144 five-meter-long 
fused silica quartz bars; the light is transmitted to the ends of the bars by total internal reflection, 
preserving the information on the angle of the light emission with respect to the track direction. 
Electrons are detected using an electromagnetic calorimeter comprising 6580 thallium-doped Csl 
crystals. Muons are detected in the Instrumented Flux Return (IFR), in which resistive plate 
chambers (RPCs) are interleaved with the iron plates of the flux return. 

The data used in the analysis were collected with the BABAR detector at the PEP-II storage 
ring. We analyze a subsample of this data consisting of 3.2 fb~^ taken on the T{4S) resonance 
and 1.2 fb~^ at a slightly lower energy to obtain a pure continuum sample. For the purpose of 
optimizing the event selection, the following simulated event samples were used: 

• 4.2 million generic T{4S) BE events generated using the Geant321 package, 

• 4.2 million continuum cc events (Geant32l), 

• 4.8 million continuum uu, dd, and ss events (Geaiit32l), 

• 50fb~^ of a fast, parametrized Monte Carlo simulation with the correct mixture of BB and 
continuum events, 

• 5,000 events for each signal mode (Geant32l). 

3 Kinematic properties of signal events 

The properties of T{4S) BB events containing a B ^ K^*^i'^£^ decay are modeled with a full 
Geant321 Monte Carlo simulation using the event generator EvtGen |^. We have implemented 
a generator based on the matrix elements computed in Ref. |jl|. Although there are significant 
uncertainties due to strong interaction effects, the matrix elements can be rigorously parametrized 
in terms of form factors, which are Lorentz-invariant functions of = M^+^_ . 

In the regions of near the J/tp and 'ip{2S) resonances, interference occurs between the penguin 
amplitude and that for charmonium production, since the J/-0 and ip{2S) can decay into . This 
interference enhances the rate on the low-mass side of the resonance and suppresses it on the high 
side. These effects are modeled in Ref. but not in our Monte Carlo generator. However, since 

are major backgrounds in our analysis, we eliminate most 
of these mass regions, as described in the following section. Small residual interference effects 
are expected for masses close to our exclusion regions, but we have ignored these in the present 
study. Further enhancements in the rate, along with related interference effects, are expected from 
charmonium resonances with masses above the tp(2S). These effects are relatively small compared 
to the expected electroweak penguin contribution, and they are ignored in this study. 

The Dalitz plot and q^ distributions resulting from these matrix elements are shown in Fig. |l|. 
For Dalitz plot variables we have chosen q^ and Ei^p, where the lepton energy is given in the B meson 
rest frame, and the positive lepton is selected for the case B^ — > K^£~^£~ and B^ — > K*£^i~ (the 
negative lepton is selected for the charge-conjugate states). For B K£~^i~ , the distribution 
peaks at low values, corresponding to rapid K recoil in the B rest frame. For B K*i~^i~ the q'^ 
distribution peaks at large values, corresponding to slow K* recoil in the B rest frame, but there 
is a pole in the amplitude at = 0, where the photon becomes real. The decay B K£~^i~ does 
not have a pole at = because B is forbidden by angular momentum conservation. 
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Figure 1: Monte Carlo simulated distributions for signal processes (a) Dalitz plot for B Ke^e , 
(b) Dalitz plot for B K*e+e-, (c) = M^+^_ distribution for B Ke+e', and (d) 
distribution for B — > K*e'^e~ . 



4 Overview of the analysis method 

The analysis method exploits the strong signature and kinematic constraints for the B K^*'>1^(.~ 
decay and the broad capabilities of the BABAR detector. Lepton and kaon identification, good 
momentum measurement, and the constraint associated with the known beam energies are essential 
for the search. Because even low-level backgrounds can be problematic, we have to be concerned 
about a possible "cocktail" of small contributions. Some of the most important sources are: 

• B ^ J/iIjK^*^ or B ^ ip{2S)K^*^ with J/i]j or 'il){2S) If one of the leptons radiates 

a photon, the mass of the dilepton system can fall below the resonance veto region. These 
events must be rejected with high efficiency, since they have precisely the same topology as 
signal events. 
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• Combinatorial backgrounds from random leptons and kaons in BB and continuum processes 
(e+e~ — > qq, where q = u, d, s, or c). In general, the lepton candidates arise from semileptonic 
decays, although leptons from photon conversions or from hadrons faking a lepton signature 
also contribute at a low level. 

• Certain backgrounds can peak in the signal region, but these are quite rare and easily rejected. 
An example is B+ D°7r+ with K+tt'. If both pions are misidentified as muons, 
this process can be confused with the signal. This type of three-body background is easily 
rejected with suitable reassignment of particle identification hypotheses and a veto on the D 
mass region. 

Our goal in the present search is to gain an understanding of background levels and distributions 
and to set limits. For this purpose, we assume that each event in the signal region is potentially 
due to the signal process and do not perform a background subtraction. However, we do carry 
through a background estimation procedure using the sidebands in data and compare the results 
with those obtained by applying the procedure to the full Geant321 simulation. 

The signal region is defined as a rectangle in the plane defined by two kinematic variables, the 
beam-energy substituted mass of the B candidate rriES) and the energy difference AE where 

WW)' 

I V^m2 +|p*|2-^/2. 
1 

The index a is over the particles that make up the candidate B meson system, nia are the masses 
of the particles and p* are their momenta measured in the T{4S) center of mass frame. \/s/2 is 
one half of the center of mass energy. It is important to distinguish between the laboratory frame, 
in which the beams have energies E^,- ~ 9.0 GeV and Eg+ ~ 3.1 GeV, and the center of mass frame, 
in which the beam energies are equal and the T(45) is produced at rest. In a symmetric-energy 
machine, these frames are identical, but at PEP-II the center of mass frame is moving with respect 
to the laboratory with a boost factor ^ 0.56. 



rriES 



\ 



AE = J 



5 Event selection 

We select events with at least five good quality tracks, with 20 or more hits in the drift cham- 
ber and originating from the beam spot. Initially, we also require either two loosely identified, 
oppositely-charged electrons, muons, or an electron-muon pair (for combinatoric background stud- 
ies). In addition, we require R2 < 0.8, where R2 is the ratio of the second and zeroth Fox- Wolfram 
moments This requirement provides a first suppression of continuum events, which have a more 
collimated ("jet-like") event topology than BB events. At this stage, R2 is evaluated in the center of 
mass frame using charged tracks only. Events passing these loose event selection criteria are further 
selected by requiring the electron and muon momenta in the lab frame to satisfy pe > 0.5 GeV/c 
and > 1.0 GeV/c. Electrons are identified on the basis of several quantities, primarily E/pe, 
where E is the energy measured in the Csl electromagnetic calorimeter and pe is the momentum 
determined by the tracking system. This ratio is near one for electrons or positrons. Muons are 
identified primarily by the number of interaction lengths of iron penetrated by the charged track 
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through the magnet flux return, which is instrumented with layers of resistive-plate chambers. 
The events are then required to lie within a large, rectangular region in the /S.E vs. rriES plane: 
mES> 5.0 GeV/c^and |AS| < 0.8 GeV. 

After these basic event selection criteria, we apply a tighter set of particle identification re- 
quirements. In addition, electrons and positrons are required to pass the conversions veto, which 
suppresses the leptons that are likely to have come from photon conversions in the material. The 
kaon identification is based on combining the energy loss information (cLE/dx) from the silicon 
vertex detector and the drift chamber for momenta p < 0.6 GeV/c, with the knowledge of the 
Cherenkov angle from the DIRC for p > 0.6 GeV/c. For the B — > Ki^i~ modes we require that 
the kaon have px > 0.6 GeV/c in the lab frame. 

For the K*^£'^i^ channels, we reconstruct the K*^ in the K^ir^ final state. The kaon 

candidate is required to be identified as a kaon, while there are no particle identification require- 
ments on the pion candidate. The mass of the K~^7r~ pair is required to be within 75 MeV/c^ of 
the K*^ mass. 

The next step is to veto B J/tltK^*^ and B tp(2S)K^*\ The requirements used for this 
purpose are somewhat complicated and are shown in Fig. |2|. We remove events with dilepton 
masses consistent with those of the J/ip or ■ip{2S). These veto regions are shown as pairs of vertical 
lines in the AE vs. M£+£- plane. However, bremsstrahlung or track mismeasurement can result 
in a large departure of the dilepton mass from the J/^l' or ip{2S) masses. Such departures are 
also accompanied by a shift in AE. We, therefore, remove these events by applying a correlated 
selection in this plane. Finally, a pair of horizontal lines is shown. There are veto regions above and 
below these lines, while the signal region lies between them. This region between the horizontal 
lines corresponds to the requirement that AE be consistent with zero. This requirement is not part 
of the charmonium veto per se but is still useful to see on this figure, because it further restricts the 
allowed region in this plane. Veto regions for the dimuon final state are defined in a similar way, 
but the definitions are somewhat different, since there is less bremsstrahlung. The B ^ J/ip K 
modes can also pass this veto if the kaon is misidentified as a lepton (most often a muon). In a 
similar way B — > D^tt, where — > K^tt^, can pass our selection criteria if both of the leptons 
are fake. These effects can be suppressed by re-assigning the particle masses and excluding mass 
combinations around the J/t/j and the D^. 

Continuum background is suppressed by using a four-variable Fisher discriminant. The variables 
are the R2, which at this stage is computed using both charged and neutral candidates, the cosine 
of the angle between the B candidate and the z axis in the center of mass frame, the cosine of the 
angle between the thrust axis of the signal and that of the rest of the event in the center of mass 
frame, and the invariant mass of the kaon and the oppositely charged lepton. The last variable 
helps discriminate against the D semileptonic decays, as it will tend to be below the D mass. Both 
continuum and BB combinatorial background is further suppressed by requiring that the combined 
B vertex probability be greater than 0.1%. 

Finally, we select the signal box in the AE vs. mEs plane: 5.272 < tues < 5.286 GeV {3a) 
and -0.10 < AE < 0.06 GeV {\AE\ < 0.06 GeV) for the electron (muon) channels. 

The net effect of our analysis selection criteria on the Dalitz plot is shown in Fig. ^. The 
requirements to remove B — > J/ipK and B ip{2S)K decays are evident as horizontal gaps, while 
those on the minimum lepton energy in the lab frame result in a loss of acceptance around the 
edges of the Dalitz plot. 
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01234 01234 

m^^^. (GeV/c^) m^^^. (GeV/c^) 

Figure 2: Regions in the AE vs. plane used to exclude J/V' K^*^ and ip{2S) K^*^ events: (a) 

shows the veto regions for the J/V'K^*) and i){2S)K^*^ events, where the J/ip or i){2S) e+e", (b) 
shows the veto regions for the J/ipK^*^ and tl'(2S)K^*^ events, where the J/ip or ip{2S) fi^ . 
The regions between the close pairs of vertical lines correspond to the nominal J/^ and ip{2S) 
resonance regions and are vetoed. The diagonal lines veto the events with bremsstrahlung and 
track mismeasurement. In the electron channel, the effects of bremsstrahlung are very substantial, 
which necessitates the use of triangular selection criteria. The signal region lies between the two 
horizontal lines except for the area excluded by the other vetos. 



6 Physics results 

Figure § shows a large AE vs. mEs region (the "grand sideband") and a small box indicating the 
signal region for each of the four modes. These scatter plots show a trend of decreasing background 
as AE increases from negative to positive values; this is a characteristic of BB background, where 
the energy of three or four randomly selected tracks is not as high as the nominal B energy in the 
center of mass frame. Events from the continuum have a more uniform distribution in AE and are 
more important for positive AE values. 

Figure |5| shows the data distribution in ruES after applying all other event selection criteria 
(including the requirement that the events lie within the AE signal region). The distributions are 
fit to the ARGUS function whose shape is derived from the large statistics fast Monte Carlo 
simulation, and are used to estimate the background contribution to the tues signal region. For 
this study, however, we do not subtract this estimated background from the yield in the signal 
region; for the purpose of setting the limits we assume that all events in the signal region might be 
due to the signal events. 

Table |l| lists the signal efficiencies, total yield, the expected background, and the 90% C.L. 
upper limits on the branching fractions. The signal efficiencies were determined from the signal 
Monte Carlo events. The efficiencies include the branching fractions for the K*^ modes. The table 
also lists the total systematic error, determined as described in Section 0. It is important to note 
that because we decided before the measurement that our result would be expressed as an upper 
limit, we do not need to apply the Feldman-Cousins procedure (Ref. fl^), for a two-sided confidence 
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Figure 3: Effects of the charmonium veto requirements on the dilepton mass distributions for (a) 



B- 



K+e+e', (b) B+ 



(c) K*^e+e-, and (d) B^ K*'^fi+fi~ 



interval, but instead use the original procedure specified in the Particle Data Book 
zero observed events corresponds to an upper limit of 2.3 events. 



in which 



7 Systematic studies 

As we are setting a conservative limit by assuming that all events observed in the signal region 
are in fact signal, the only systematic effects that affect the limits are those due to the modeling 
of our signal efficiency and the number of produced B mesons in our sample. Table |2| summarizes 
the systematic uncertainties that we have considered. The uncertainty on the tracking efficiency 
is 2.5% per track. Lepton and kaon identification efficiencies are measured in control samples in 
data (e.g., a Bhabha sample for electrons), where a measure of the uncertainty is obtained by 
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Figure 4: vs. m£;5 (grand sideband) for data: (a) B'^ K'^e'^e , (b) 

B^ — > K*^e^e~ , and (d) B^ K*^fi^fi~. The smaller boxes show the signal region. 



comparing different control samples. We assign an error of 2.0% per electron, 2.5% per muon, and 
3.0% per kaon. We have also studied the effects of momentum smearing and shifting on the AE 
selection criteria, and assign an error of about 3.0% for this effect. Comparing data and Monte 
Carlo efficiencies for the mES requirement in the J/ip control region, we assign a 3.0% error for 
the modeling of this requirement. We have compared the data and Monte Carlo efficiencies for the 
B vertex probability and the Fisher discriminant selection criteria. For the former, we assign a 
systematic error of 3.0% for the B~^ — > K^£'^£~ modes and 4.0% for the B^ K*^i'^£'' modes. 
For the latter, we assign an error of 3.0% for all of the four modes. Furthermore, as a check that 
we understand our signal efficiencies, we have compared the yields in data and Monte Carlo events 
for the B+ J/^K+ and B^ J/ijjK*^ (with ^ K+tt") channels. We have used the world 
average Q branching fractions for B^ — > J /ipK'^ and B^ J /ipK*^ and found agreement within 
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Figure 5: mES for data after all other event selection criteria are applied: (a) K'^e^e , (b) 

B+ K+n+fi-, (c) B° K*^e+e-, and (d) K*^ . The shape of the fit (the ARGUS 

function) is obtained from the large statistics sample of fast parametrized Monte Carlo events. The 
lines indicate the signal region. 



the statistical errors between the data and Monte Carlo samples. 



8 Summary 

We have searched for rare B decays B K^*^ £+£- in a sample of 3.7 x 10^ BB events. We find no 
observable signal for any of the four modes considered, and set preliminary 90% C.L. upper limits 
on the branching fractions of 

B{B+ K+e+e-) < 12.5 x 10"^ 
B{B+ K+fi+fi-) < 8.3 X 10"^ 
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Table 1: Signal efficiencies, systematic uncertainties (combining the uncertainties on the signal 
efficiencies and on the number of produced T(4S) mesons, as listed in Table ||), the number of 
observed events, the number of estimated background events, and upper limits on the branching 
fractions. In computing the upper limits we have assumed B{K*^ K+n-) = 2/3. 



Mode 


Efficiency (%) 


Total systematic 


Observed 


Bkgd. estimated 








uncertainty (%) 


events 


from data 


(90% C.L.) 


B+ K+e+e- 


13.1 


11.7 


2 


0.20 


< 12.5 


B+ K+n+fi- 


8.6 


12.3 





0.25 


< 8.3 


_^ K*^e+e- 


7.7 


14.2 


1 


0.50 


< 24.1 




4.5 


14.8 





0.33 


< 24.5 



Table 2: Summary of the systematic uncertainties on the signal efficiencies and the number of 
produced T(4S) mesons as a percentage error on the branching fraction, B. The total systematic 
error is the sum of the individual contributions added in quadrature. 





(AB/B) (%) 


Kee 




K*ee 




Tracking efficiency 


7.5 


7.5 


10.0 


10.0 


Lepton identification 


4.0 


5.0 


4.0 


5.0 


Kaon identification 


3.0 


3.0 


3.0 


3.0 


AE requirement efficiency 


2.0 


2.5 


3.3 


1.5 


mES requirement efficiency 


3.0 


3.0 


3.0 


3.0 


Vertex requirement efficiency 


3.0 


3.0 


4.0 


4.0 


Fisher requirement efficiency 


3.0 


3.0 


3.0 


3.0 


Number of produced T(4S) 


3.6 


3.6 


3.6 


3.6 


MC signal statistics 


3.4 


4.0 


4.3 


6.1 


Total systematic error 


11.7 


12.3 


14.2 


14.8 



B{B^ K*°e^e-) < 24.1 x 10"^ 
B(S° ^ K*^n+fi-) < 24.5 X 10~^ 

The limits for the B~^ K'^l'^l^ modes are comparable to those set by other experiments, while 
those for B^ K*^£'^i~ are less sensitive with this data sample. We plan to analyze substantially 
more data in the near future. 
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